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We report our recent studies on the finite density QCD obtained from lattice QCD simulation with clover- 
improved Wilson fermions of two flavor and RG-improved gauge action. We approach the subject from two paths, 
i.e., the imaginary and real chemical potentials. 

1 Introduction 

QCD at finite temperature and density has been one of the most attracting subjects in physics. Many phenomeno- 
logical models predict that the QCD phase diagram has a very rich structure, and thoroughgoing analyses of heavy 
ion data show that we are sweeping finite temperature and density regions. See Ref. HO. 

First-principle calculations based on QCD are now highly called. If such calculations would be at our hand, 
their outcomes are also very valuable for many research fields: high energy heavy ion collisions, the high density 
interior of neutron stars and the last stages of the star evolution. Needless to say, the inside of nucleus is also a 
baryon rich environment, and lots of contributions to nuclear physics could be expected. 

Unfortunately, the first principle lattice QCD simulation suffers from the sign problem. Nevertheless, there 
have been many progresses such as the reweighting method, the imaginary chemical potential and the canonical 
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formulation; now some light is shed on the QCD phase diagram. For reviews, see e.g. [2, 3]. 

Here, we report our recent trials to promote the finite density lattice QCD. It contains two results (4. 5\: (he 
determination of the phase boundary of the deconfinement transition based on the imaginary chemical potential 
approach and a reduction formula for the Wilson fermion determinant. 



2 Imaginary Chemical Potential Approach 
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Figure 1: Schematic figures for the Nf = 2 QCD phase diagram in the (/.t 2 , T) plane (left) and (/xj/T, T) plane 
(right). A : Pseudo-critical point at \i = 0. B : Critical endpoint. C : Roberge-Weiss endpoint. AB : Pseudo-critical 
line. AC : Extension of the line AB into the imaginary chemical potential plane. CD : Roberge-Weiss phase 
transition line /Zj/T = tt/3. In the right panel, larger /i/ /T region of the phase diagram is obtained from the RW 
periodicity. 



The QCD with an imaginary chemical potential is free from the sign problem. Using a relation 

(det A(/x))* = det A(/i), (fi = fi R + ifn), 



(1) 



it is straightforward to prove that det A(/i) is real for fi = ifij. A partition function and its free-energy are analytic 
within one phase even if chemical potential is extended to complex, which is true until the occurance of a phase 
transition. This validates the imaginary chemical potential approach for the study of the QCD phase diagram. 
In addition, the QCD phase diagram in the imaginary chemical potential regions have a unique feature called 
the Roberge-Weiss periodicity J3l, see Fig. [T] There have been several studies in staggered fermions I7 - J_6 ] and 
standard Wilson fermions 11711 . 
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We employ a clover-improved Wilson fermion action of two-flavors and a renormalization-group improved 
gauge action. The clover-improved Wilson fermion action is given by 

A(i,») = 6 X , X , -kJ2 [(1 -li)Ui{x)8 x , iX+l + (l+7^(z')<Ws" 

i=l 

-« [e +Al (l - 74)^4(2;)^, x+l + e-"(l + 74)C/ 4 t (a ; ')^',x-4] " ^ c swS x , x > ^ 

Here fi is the quark chemical potential in lattice unit, which is introduced to the temporal part of link variables. 

In order to scan the phase diagram, simulations were done for more than 150 points on the (/xj, f3) plane in the 
domain < fii < 0.28800 and 1.79 < j3 < 2.0. All the simulations were performed on a 7V S 3 x N t — 8 3 x 4 lattice. 
The RW phase transition line in the present setup is given by [ii = tt/12 ~ 0.2618. The value of the hopping 
parameter k were determined for each value of f3 according to a line of the constant physics with nips / ray = 0.8 



obtained in Ref. 118 1. 

Scatter plots of the Polyakov loop in the complex plane are shown in Fig. [2] where we choose two typical 
cases ft = 1.80 for the hadronic phase and f3 — 1.95 for the QGP phase. At low temperatures, the Polyakov 
loop is small in magnitude for any /i/ and continuously changes in a clockwise direction as increasing [ij. On the 
other hand, at high temperatures, the Polyakov loop grows to 0.2 ~ 0.3. It stays at the real axis for [ij < n/12 
and jumps to the left-lower side at [ij = 7r/12. The difference of the Polyakov loop modulus between high and 
low temperatures shows the deconfinement crossover, which is the curve AC in Fig.Q] The observed jump of the 
Polyakov loop at fij = 7r/12 is the Roberge- Weiss phase transition, which is the line CD. Thus, the phase structure 
in fj, 2 < regions of the QCD phase diagram can be determined by observing the behavior of the Polyakov loop. 
The properties of the phase transitions and RW endpoint are obtained from the susceptibility of the Polyakov loop. 
We obtain the location of the RW endpoint (3 = 1.927(5), which corresponds to T/T pc ~ 1.15. 

Critical values of (3 for the deconfinement transition are obtained from the susceptibility of the Polyakov loop 
modulus for each /i/. Using the data for the critical values of j3, we can determine the pseudo-critical line. Obtained 
pseudo-critical line is analytically continued to fi 2 > region. Th results are shown in Fig. [3] where we employ 
physical unit (fi = fia). The curvature at jX/T pc — of a power series of {jl/nT pc ) 2 is £2 = K 2 d,2 = 0.38(12). 



The present results are slightly smaller than other studies, see e.g. Ref. (119| 

3 Reduction Formula for Wilson Fermions 

In the lattice QCD simulations with finite chemical potential /j, often we must handle the fermion determinant 
det A(/i), directly. For example, the reweighting method requires a ratio of two determinants, d°t A*(>) ' ^ e ^ en " 



sity of state method needs the phase information [20]; The canonical formulation needs the Fourier transformation 
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Figure 2: Scatter plots of the Polyakov loop. Left : ft — 1.80 (low temperature (below T pc )). Right : j3 = 1.95 
(high temperature (above Trw))- 




Figure 3: The pseudo-critical line j3 pc in the imaginary (left panel) and real (right panel) region. 
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of the fermion determinant. In these approaches, the heaviest part of the numerical calculations is the evaluation 
of the determinant. An efficient way of the determinant evaluation is highly desirable. Here we introduce a matrix 

n 

reduction formula for Wilson fermions, which was first constructed by Borici [21]. Later it was studied with the 
inclusion of the fugacity expansion 22L 

The Wilson fermion matrix defined in Eq. (ffj can be divided into three terms according to their time dependence 



A = B - 2z~ 1 Kr-V - 2znr + VK 



Here r± = (r ± 74)72 with the Wilson parameter r and z = e p . Each component is defined by 



(2) 



B 



o 

(*,*') = 5 X , X , -n^2{(r - li)Ui{x)8 x , tX+ . + (r + ^)Uj(x')6 x ,^ + S C u 



V(x,x') = U 4 (x)S x , x+h V*(x,x') = f/i(^)^,x-a- 



(3) 



(4) 



They satisfy VV^ = I. Note that r± are projection operators in the case that r = 1. 

Now, we construct a reduction formula for the Wilson fermions. A starting point is to define a permutation 
matrix P = (c a r_ +cj 3 r + Vz~ 1 ) J21]]. The parameters c a and q, are arbitrary scalar except for zero, and may be set 
to one. Since r± are singular, the matrix P must contain both of them; otherwise P is singular. It is straightforward 
to check det(P) = {c a CbZ~ x ) N I 2 , where N = AN c N x N v N z N t . Multiplied by P, the quark matrix is transformed 
into 



AP = (c a Br_ - 2c b nr + ) + (c b Br + - 2c a nr_)Vz~ 1 . 



(5) 



Carrying out the temporal part of the determinant, we obtain 



det AP 



fhz- 1 
a 2 faz^ 1 



«3 



= ^JJ det( ai )^ det (1 + z~ Nt Q) , (6) 
where Q = (a^ 1 f3\) ■ ■ ■ (a^/3/v t ), which is often referred to as a reduced matrix or transfer matrix. The block- 



5 



matrices a and j3 are given by 

a i = a ab '^(x,y ) t i ) 

= c a B ab ^(x, y, U) r<L v - 2c b n r^5 ab S(x - y), (7) 
(3 t =p ab ^(x,y,U), 

= c b B ac ^{x, y, ti) r?U?(y, i,) - 2c a n r^6(x - y)U?{y, U), (8) 

where the dimensions of a.j and are given by iV rc( j = N/N t — 4N x N y N z N c . Substituting det(P) = 
(c a CbZ~ 1 ) N l 2 , we obtain 

det A = {c a c b )- N ' 2 z- N ' 2 det ^JJ a)j det (z w * + Q) . (9) 

Here, Q is independent of /i and its rank is given by iV re d = N/N t , while that of the Wilson fermion is originally 
given by N. 

With the eigenvalues A„ = {A| det(Q — XI) = 0}, the determinant of Q is given by 

det{z N * +Q) = Y[(\ n + z Nt ). (10) 
n=l 

Expanding this in powers of the fugacity finally obtain the reduced quark determinant 

A f rod/2 

detA(/x)= C4e M/T )", (11) 

n=-Af rod /2 

Note that we redefine the index of c n (c n by c_ ra ) to obtain the second line from the first one. Here, C n = Cc n 
with C = (c a c b )- N / 2 (nfi det(aO). 

Using a relation Eq. (Q~|l and the reduction formula, one gets 

N Icd JV rcd 

71 — 1 71 — 1 

where £ = 2r *. This holds for any £ G C. For £ = — A„, the left-hand side vanishes, and so should be the 
right-hand side. Then, the eigenvalue always appear in a set 

An, I/A; (13) 



The relation is also pointed out by Alexandru and Wenger II22I1 . 
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Figure 4: Schematic figures for the reduction procedure. 




Figure 5: Distributions of the large eigenvalues in complex plane. 
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Figure 6: Distributions of the small eigenvalues in complex plane. 
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Figure 7: Histogram of the absolute value of A. Left : large eigenvalues. Right : Small eigenvalue. 
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Figure 8: Histogram of the phase of A. 
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The reduction formula makes it easier to calculate fermion determinant. We plan to evaluate the phase transition 
line at real chemical potential points. Combining estimations of the phase transition line both at real and chemical 
potential regions, we will get more reliable information about QCD phase structure. 

Figures [5] and [6] show the scatter plot of {A} for three different temperatures /3 = 1.80,1.855,2.0, which 
correspond to T/T c — 0.9, 1, 1.3, respectively. Figures|7]and[8]show the histogram of the eigenvalue distribution. 
The simulation setup was the same as given in the previous section. Note that Fig. [6] enlarges a small domain near 
the origin in Fig. [5] 

The eigenvalues are distributed in two separate regions, and there is a margin between the two regions, where 
no eigenvalue is found. The histogram of the absolute value of A also show this behavior, see the right panel of 
Fig. El 

The /3 dependence appears in the phase of A; A are distributed in a Z% symmetric manner at low temperatures, 
while not at high temperatures. The symmetric property is broken at high temperatures, and A approach to real 
axis. Note that this behavior is observed both for small and large eigenvalues, because of the pair nature of A. 

The properties of the fugacity coefficients c„ was discussed in |4j]. 

The simulation was performed on NEC SX-8R at RCNP, and NEC SX-9 at CMC, Osaka University, and 
HITACHI SRI 1000 and IBM Blue Gene/L at KEK. This work was supported by Grants-in-Aid for Scientific 
Research 20340055 and 20105003. 
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